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r^ ^ Abstract. Scale-invariant fluctuations of growing interfaces are studied for 

^ , circular clusters of an off-lattice variant of the Eden model, which belongs to 

^ ■ the (1 -I- l)-dimensional Kardar-Parisi-Zhang (KPZ) universality class. Statistical 

a' properties of the height (radius) fluctuations are numerically determined and 

' compared with the recent theoretical developments as well as the author's 

experimental result on growing interfaces in turbulent liquid crystal [K. A. 
Cw , Takeuchi and M. Sano, arXiv:1203.2530]. We focus in particular on analytically 

TT ■ unsolved properties such as the temporal correlation function and the persistence 

probability in space and time. Good agreement with the experiment is found 
|— ^ ' in characteristic quantities for them, which implies that the geometry-dependent 

e universality of the KPZ class holds here as well, but otherwise a few dissimilarities 

are also found. Finite-time corrections in the cumulants of the distribution are 
I , also studied and shown to decay, for the mean, as i~^" within the time window 

^ 1 of the simulations, instead of t~^i^ which arises as the typical leading term in the 

Ch ' previously known cases. 

! PACS numbers: 05.70.Jk. 02.50.-r, 89. 75. Da, 81.15.Aa 

> 

CO ' 1. Introduction 

00 , 

\l , Growing interfaces driven by local stochastic interactions constitute a prototypical 

example where scale invariance leads to universality even in systems far from 



(N 



CO . equilibrium [1]. In their seminal work [2], Kardar, Parisi, and Zhang (KPZ) proposed 



o 

(N 



a continuum equation for describing such growing interfaces, which reads 

^^h{x,t)^v^^h+^{^hf+i{x,t) (1) 

, , , with the local height h{x^ t) of the interfaces and the white Gaussian noise {£,{x, t)) = 

rS ■ and (^(a;,t)^(a;',i')) = D5{x — x')S{t — t'). They derived a set of characteristic 

jrt ! exponents a = 1/2, /3 = 1/3 and z = a/(3 — 3/2 that describe scale-invariant 

fluctuations of the interfaces in 1 -I- 1 dimensions, in agreement with the earlier 
theoretical studies on related problems [3,4]. These exponents are universal as widely 
confirmed by theoretical and numerical work [1] as well as by a few experiments (see [5- 
7] and references therein), constituting the KPZ universality class. Recently, studies in 
this context entered an unprecedented phase, when Johansson rigorously derived the 
asymptotic distribution function for a model in the KPZ class [8] on tha basis of related 
combinatorial problems [9]; this spearheaded remarkable theoretical achievements 
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marked in the last decade (for recent reviews, see [10-12]). Their main conclusions 
are twofold: (i) The distribution function and the spatial correlation function were 
obtained analytically for a few solvable models, including the KPZ equation [13-19], 
and revealed nontrivial connection to random matrix theory, (ii) These are expected 
to be universal, but nevertheless depend on the global shape of the interfaces, e.g., 
whether they are curved or flat, as first pointed out by Prahofer and Spohn [20]. 
In particular, the asymptotic distribution function for the curved interfaces is given 
by the largest-eigenvalue distribution of large random matrices in Gaussian unitary 
ensemble (GUE) [8,20], called the GUE Tracy- Widom distribution [21], and for the 
flat interfaces by the corresponding distribution for Gaussian orthogonal ensemble 
(GOE) [20]. 

This detailed yet geometry-dependent universality of the KPZ class was recently 
underpinned by an experiment on growing interfaces of topological-defect turbulence in 
nematic liquid crystal, in which the author and his coworker identified the analytically 
obtained distribution functions and the spatial correlation functions for both circular 
and flat interfaces [5-7] . In [7] they further advanced their analyses to study statistical 
properties that remain analytically unsolved even for the solvable models, in particular 
those related to the temporal correlation of the interface fluctuations. This revealed 
that the difference between the circular and flat interfaces is not restricted to the 
distribution and spatial correlation functions, but actually arises in many other 
quantities as well, sometimes even with qualitative differences, as summarized in table 
1 of [7]. These quantities — specifically, the temporal correlation function and the 
temporal and spatial persistence probabilities — are expected to be universal in the 
asymptotic limit when appropriately rescaled. 

Here, as an approach complementary to the experiment, the author performs the 
same analyses for a numerical model, namely an ofF-lattice variant of the Eden model 
[1,22] introduced in the present paper, and provides a numerical case study on the 
universality in those unsolved quantities. Focus in the present study is set on the 
circular interfaces, which are numerically much less studied than the flat case and are 
also more delicate: because the circumference, or the system size of the interface, is by 
construction finite at finite times and grows with time, the ensemble and the spatial 
averages lead to different results when they are used to define the statistical quantities 
[7] . Using the spatial average instead of the ensemble one yields false estimates even 
for the growth exponent /3 [7]. Since this point has not been explicitly recognized in 
past studies, it is important to provide a set of reliable numerical estimates for the 
circular interfaces of a numerical model — this is the aim of the present paper. 

2. Model 

Only few numerical models are known to produce circular interfaces, most of which 
are variants of the Eden model. The Eden model was originally defined on a lattice. 
Starting with a seed particle at the origin, at each time step one adds a new particle on 
a randomly chosen site on the perimeter of the cluster [1, 22]. This however turned out 
to induce anisotropy due to the lattice structure [1, 23] unless a specific growth rule is 
introduced, such as the growth probability dependent on the number of the occupied 
nearest-neighbors [24,25]. Off- lattice versions of the Eden model have therefore been 
considered occasionally in the literature [25-28] , but in most cases time t is measured 
by the global radius {h) of the radial cluster, which is not correct at finite times 
because the two quantities are in general connected by d{h)/At ~ t^oo + c^t'"^^"^ [5- 
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Figure 1. Off-lattice Eden model, (a) Sketch of the model. A new particle 
(dashed) is being added to its ancestor (shaded). In Alves et al.'s model [25,27, 
30], a random number for the angular position of the new particle is generated 
within [tt/G, — vr/6] , where (p is defined in the sketch and 7r/6 is due to the 
finite radius of the new particle. In contrast, in our model, the random number is 
always drawn from [0, 2it) and then the attempt of generating the new particle is 
accepted with probability {{J> — tt/3)/2it. (b) Typical cluster produced by our Eden 
model. The active particles, or the interface, at t = 405, 803, 1174, 1592, 1999 are 
shown. 



7,29]. To resolve these problems, as well as to study a numerical model that the 
author considers better corresponds to a coarse-grained picture for the experimentally 
investigated growth of the liquid-crystal turbulence [5-7], a new variant of the ofF- 
lattice Eden model is introduced and studied in the present paper. 

The model is defined as follows. First, place a round particle of unit diameter at 
the origin of two-dimensional continuous space. Identical particles are then added one 
by one according to the following steps, where N^ is the number of the active particles 
as defined below. The initial particle is active, hence A^a = 1 at i = 0. (1) Choose 
one of the N^ active particles randomly. (2) Attempt to place a new particle on its 
border, in a direction randomly chosen in the range [0,2tt). (3) Place this particle if 
it does not overlap any other particles, otherwise abandon the attempt [figure 1(a)]. 
(4) Label as inactive those particles to which no particle can be added any more for 
lack of empty space, thereby excluding them from the list of the N^, active particles. 
Since we are interested in the dynamics of the interface, we also exclude the particles 
left in the bulk, which are surrounded by an outer closed loop of adjacent particles. 
Here, particles are regarded as adjacent if the distance between them is shorter than 
V3 because in this case an inner particle cannot infiuence outer ones in any manner. 
Note also that the inactive particles remain obstacles for the active ones; in practice, 
this is realized by recording the angles of the inactive neighborhood for each active 
particle. (5) Increase time t by 1/iVa, whether the new particle is placed or not. This 
choice indicates that each particle has a chance to add a neighbor once per unit time 
on average. Therefore, it is statistically equivalent to use the whole particle number 
A'' instead of A"a for all the steps above, but the use of the active particle number A'a 
speeds up the simulations significantly. 

This model resembles the off-lattice Eden model studied by Alves et al. [25] (also 
described in [27,30]) except for the criterion of the inner particles and for the way a 
new particle is added to the randomly chosen ancestor. On the former point, Alves et 
al. inactivated the particles inside a central core of radius 0.8f [27, 30], where f is the 
mean distance of the particles on the interface from the origin, whereas in our model 
it is determined by the purely geometrical consideration, for which there is no risk 
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of false inactivation. The latter point concerns the step (2) in the above procedure. 
When determining the direction to place a new particle, Alves et al. drew a random 
number from such angles that correspond to the empty neighborhood of the ancestor, 
while we always choose an angle within [0, 27r) and then judge whether the particle can 
be placed or not [figure 1(a)]. This guarantees local isotropy of the growth process. 
The author also considers that, in coarse-grained scales, it better corresponds to the 
growth of topological defects in the liquid-crystal turbulence [5-7] . 

The data presented below are obtained from 3000 independent simulations of the 
above-defined Eden model up to i = 1999, which roughly corresponds to the final 
cluster size of iV = 2 x 10^ if one does not discard the inner particles. Figure 1(b) 
shows a typical growing cluster; it displays the active particles at different times, 
labeled by i below, which form the interface at each time by construction. The local 
interface height /i(a;, t) is then determined from the distance r^ between the origin and 
the active particle i at a given angular position. More specifically, the full azimuthal 
range [0, 2tt) is divided into [27r(ri)J bins, where (r^) is the ensemble average of r^, or 
the mean radius of the clusters at the given time, and [• • -J represents the integer part; 
then, for each bin and cluster, a single value of h{x, t) is assigned as the mean of ri 
within the corresponding range of angular positions. Note that the width of the single 
bin is then roughly equal to the diameter of the particles. The lateral coordinate x is 
defined along the mean shape of the clusters, i.e., the circle of circumference 27r(ri). 

3. Results 

We first check the scaling exponents a and /3 by the usual method based on the 
Family- Vicsek scaling [31], which describes the spatial and temporal dependence 
of the roughness growth. To this end, we define the interface width w{l,t) as the 
standard deviation of the height h{x,t) measured over length /, or, more specifically, 

w{l, t) = ( J {[h{x, t) — {h{x, t))i\ )i ) with the ensemble average (• • •) and the average 

(• • •); taken inside a segment of length I around the position x. The Family- Vicsek 
scaling then reads 

w{l,t)^t^FUlt''/n-\. , , , (2) 

[tp ioTi:$>h, 

with z = a//3, a scaling function F^ and a crossover length scale /, ^ t^'^. This 
relation is clearly confirmed in our numerical data, together with the KPZ-class 
exponent values a — 1/2, (3 — 1/3 and z = 3/2 (figure 2). The width w{l,t) at 
fixed times grows as w ^ /" with a = 1/2 for short length scales /, while it saturates 
at the value that increases with time [figure 2(a)]. This width at the largest length 
scale can be accurately measured by the overall width W{t) = ^J{[h{x,t) — (h)]"^), 
which indeed shows a clear power law W ^ t^ with /? = 1/3 [figure 2(b)]. Note here 
that the overall width W{t) should be measured with the ensemble average (• • •) and 
not the spatial average taken for each interface, which turned out to bias the apparent 
value of /3 as detailed in [7]. The agreement with the KPZ-class exponents is also 
checked by data collapse: the data for the width w{l,t) at different times collapse 
reasonably well onto a single curve F^ when wt~^ is plotted against lt~^'^ [insets of 
figure 2(a)]. 
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Figure 2. Family- Vicsek scaling, (a) Width w{l,t) against length I at different 
times t = 24, 82, 322, 1999 (from bottom to top). The inset shows the same data 
in the rescaled axes, wt~^ vs lt~^'^, with the KPZ-class exponents /3 = 1/3 
and 2 = 3/2. (b) Overall width W(t) against time t. The inset confirms the 
convergence to the asymptotic power law W ^ t^ with /3 = 1/3. The dashed lines 
are guides for the eyes indicating the slopes for the KPZ-class exponents. 



The realization of the KPZ-class growth exponent /? = 1/3 implies that the local 
height h of the interfaces can be described with a deterministic linear growth term 
and a stochastic t^'^ term, as follows: 



h^v^t+{Ttf'^X, 



(3) 



where tioo and T denote two constant parameters and x * random variable that 
captures the fluctuations of the growing interfaces. Note that we focus for the 
moment on asymptotic one-point statistics of the interface fluctuations and hence 
do not consider the dependence of x on a: and t. 

First we estimate the values of the two parameters Wjx) and F, in the same 
manner as for the experiment on the liquid-crystal turbulence [5-7]. The linear 
growth rate Voo is obtained as the asymptotic growth speed of the mean height. More 
specifically, since equation (3) indicates d(/i)/di ~ Vao + Cyt~'^'^ with a constant c„, 
linear regression of the data for d{h)/dt against f^^/'^ in the asymptotic regime, here 
t > 250, provides a precise estimate of Woo at Woo = 0.5139(2) [figure 3(a)]. Concerning 
the amplitude F for the t^/'^-fluctuations, we use the relation to the second-order 
cumulant: {h'^)^ = W{t)'^ ~ (ri)^^^(x^)c- Setting the arbitrary variance of x to that 
of the compared distribution function, namely the variance of the GUE Tracy- Widom 

distribution denoted by (Xgue)c' ^® P^*-** ^^ figure 3(b) ((^^)c/(Xgue)c)'^^^/^ ^^^ 
read its time asymptotic value. It shows a considerable finite-time effect, which is also 
visible in the overall width W{t) [figure 2(b)]. Following the method used for the fiat 
interfaces in the liquid-crystal experiment [7], we first fit a power law ait^ +02 to the 
data for t > 10 in figure 3(b) and obtain S = 0.70(3). Since it is closest to 2/3 among 
the multiples oi (3 = 1/3, we assume 5 = 2/3 and fit again the data with a weight 
proportional to (y — a2)^^, where y is the ordinate of the data as shown in figure 
3(b). This finally provides our estimate at F = 1.01(2). Note that the parameter 
values obtained here are connected to the three parameters in the KPZ equation (1) 
by t;o„ = A [7,29] and F = A21A1/2 with A = D/2v [10]. 
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Figure 3. Parameter estimation, (a) Estimation of the linear growth rate Vao- 
The time derivative of the mean height, d(h)/dt, is plotted against t~^'^. The 
y-intercept of the linear regression (dashed line) provides the estimate of v^c ■ (b) 
Estimation of the amplitude F for the i^' "^-fluctuations. The time-series of the 
second-order cumulant (?i^)j, is displayed in such a way that the value of F is 
found in the limit t — >■ oo. It is obtained by fitting at~^''^ + F (dashed line) to 
the data for t > 10 in order to take into account the finite-time effect (see text). 



Using the measured parameter values, we rescale the height h{x, t) in such a way 
to extract the random variable x'- 

Figure 4(a) displays the histograms of the rescaled height q at different times (symbols) 
and compare them with the GUE and GOE Tracy- Widom distributions (dashed and 
dotted lines, respectively). The data at the later time clearly indicate the GUE 
Tracy- Widom distribution, in agreement with the recent analytical results for the 
curved interfaces in the solvable models [10-12]. This is more quantitatively checked 
by plotting the difference in the nth-order cumulant, (9")^. — (Xgue)c' ^'^ figure 
4(b). We thereby confirm that the circular interfaces of our off-lattice Eden model 
indeed belong to the KPZ universality class at the level of the one-point distribution 
function. Moreover, the differences in the nth-order cumulants also measure the 
finite-time corrections toward this asymptotic distribution, which turn out to decay as 
\{'f')c ~ (Xgue)cI "" t~'^^^ for n = 1 and 2 [figure 4(c,d)]. The power-law decay i"^/^ 
for the second-order cumulant has also been reported in all the previously studied 
systems [7,13,14,32,33]. In contrast, noteworthy is the decay i~^/'^ found for the 
first-order cumulant, or the mean, because the typical leading correction shown by 
the previous studies is i"-^/"^, both for the solvable models [13, 14,32,33]| and for the 
liquid-crystal experiment [5-7]. On the one hand, one cannot exclude the possibility 
that the correction in our case shows a crossover from t~'^/^ to t~^^^ at larger times, 
as reported in simulations of flat interfaces in the discretized KPZ equation [34]. On 
the other hand, our result may indicate that our Eden model is endowed with some 

X Ferrari and Frings [32] argued that the leading correction term for the mean should be in general 
0{t~^/^) for discrete models, with mathematical demonstrations for the polynuclear growth (PNG) 
model and the totally and partially asymmetric simple exclusion processes (TASEP and PASEP). 
They, however, also showed a couple of exceptions where this term vanishes, namely the flat PNG 
interface and the PASEP when the hopping rate is tuned to the critical value [32]. 
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Figure 4. One-point distribution and finite-time corrections, (a) Histogram of 
the rescaled local height q = (h — Voot) / (Tt)^ ' ^ for t = 82 (blue circles) and 
t = 1999 (red diamonds). The dashed line indicates the GUE Tracy- Widom 
distribution expected for the curved interfaces, whereas the dotted line shows the 
appropriately rescaled GOE Tracy- Widom distribution for the flat interfaces, (b) 
Finite-time corrections in the nth-order cumulants. Shown are the differences 
between the cumulants of the measured rescaled height q and the theoretical 
values for the GUE Tracy- Widom distribution. (c,d) Corrections in the first- and 
second-order cumulants [same data as (b)] shown in the double logarithmic scales. 
The dashed lines are guides for the eyes, both indicating the slope —2/3. 



kind of symmetry, for which the leading correction term t^^'^ for the mean vanishes. 
Note also that Alves et al. [25] numerically found the GUE Tracy- Widom distribution 
in their off-lattice Eden model, as well as in two on-lattice isotropic Eden models, and 
reported finite-time corrections for the mean proportional to i^^/'^. This difference 
should arise from the different definitions of the off-lattice Eden model as described 
in the previous section. 

The spatial correlation function 



Cs{l;t) = {h{x + l,t)h{x,t)) - {h{x + l,t)){h{x,t)). 



(5) 



is also a quantity of central interest in the recent analytical studies [10-12]. For the 
solved cases of the curved interfaces, it has been shown to coincide with the covariance 
of the stochastic process called the Airy2 process A2{t') [35], which was noticed to be 
equivalent to the dynamics of the largest eigenvalue in Dyson's Brownian motion for 
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Figure 5. Spatial correlation function Cs(i;t). (a) Rescaled correlation function 
C^{C;i) = Cs(«;t)/(rt)2/3 against rescaled length C = (AZ/2)(rt)-2/3. The 
symbols show the results for the Eden model at t = 82 (blue circles) and 
t = 1999 (red diamonds), while the dashed line indicates the Airy2 covariance 
92 (C)- (b) Finite-time correction in the integral of the rescaled correlation function 
Cs"*(*) = Io° C'sCCi *)dC- Displayed is the difference from the value for the Airy2 
covariance gJ,"' = f,°° g2(C)dC as a function of time. The dashed line is a guide 
for the eyes showing the slope —1/3. 



a{l;t)^{Ttf/^gJ^{Ttr'/' 



(6) 



GUE random matrices [36]. The prediction reads 

Y 

with 52(C) = {■^2{t' + 0-^2{t')), which was indeed substantiated in the liquid-crystal 
experiment [5,7]. Our numerical data also confirm this in the time asymptotic 
limit, as shown in figure 5(a) which displays the rescaled correlation function 
C^(C;t) = Cs{l;t)/{rt)^/^ against C = {Al/2){Tt)-^/^. The finite-time correction 
is also quantified by measuring the integral Cg"'(i) = L C!,{(^;t)dC as a function of 
time. This indeed approaches the value for the Airy2 covariance g™^ = /„ g2{C)d-C 
by a power law 32 — C's"'(^) ^ t~^'^ [figure 5(b)], which was also identified in the 
liquid-crystal experiment [5,7]. 

From now on we study statistical quantities that remain out of reach of 
rigorous theoretical treatment so far. Among them particularly important are those 
characterizing the correlation in time, especially the temporal correlation function 

Ctit,to) = {h{x,t)h{x,to)) - {h{x,t)){h{x,ta)). (7) 

The temporal correlation should be measured along the directions in which fluctuations 
propagate in space-time, called the characteristic lines [37,38]. In our case these are 
simply the radial directions of the circular growth and represented by the fixed x in 
equation (7). 

For the circular case, Singha [39] performed an approximative theoretical 
calculation based on the assumption of effectively linear evolution of the height 
fluctuations, and derived 
Ct{t,to) 






Fsingha.{t/to'ib), 

_ e''(i-i/v^)r(2/3,&(l 



1/Vi)) 



r(2/3) 



(8) 
(9) 
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Figure 6. Temporal correlation function Ct(t, to)- (a) Ct{t,to)/Ct{to,to) against 
t/to with different to (as shown in the legend). The dashed lines show the results 
of fitting by Singha's function (9). Notice that the correlation function remains far 
above zero, (b) The values of the fitting parameter b in equation (8) as a function 
of ifl. The dashed line indicates the mean value of b estimated for t > 500, while 
the dotted curve shows the result of a power-law fit fe(to) = "-to 4 



displays the quality of this power-law fit in the axes 
for the eyes indicating the slope for S = 1/2. 



)oo- The inset 
6 vs 4oi with a guide 



with a single unknown parameter b, the upper incomplete Gamma function r(s,x) = 
J^ y^~^e~ydy and the Gamma function T{s) = r(s,0). For the liquid-crystal 
experiment [7], this functional form works only asymptotically. More specifically, it 
fits the experimental data for t > to on condition that the right hand side of equation 
(8) is multiplied by a time-dependent coefiicient c(to), which turned out to satisfy 
c{to) — )■ 1 with to ^ CO and is attributed, somewhat speculatively, to microscopic 
dynamics of the liquid-crystal convection. On the other hand, Singha himself checked 
his prediction by simulations of the conventional on-lattice Eden model [39] . He then 
found good agreement if an independent value of b is chosen for each growth direction, 
without however studying the dependence on io- 

The temporal correlation function is measured in our off-lattice Eden model and 
shown in figure 6(a), together with the results of fitting by Singha's function (9). As 
shown in the figure, Singha's function fits our data reasonably well within the whole 
time window, or c(io) = 1 for all to-, in agreement with his simulations for the on-lattice 
Eden model. We however find that the value of the parameter b obtained by the best 
fitting of the data increases with to [figure 6(b)]. Assuming the convergence of b at 
later times such as i > 500, we estimate the asymptotic value of 6 at & = 1.18(3). 
On the other hand, attempt of a power-law fit in the form b{to) ~ at^ + b^o yields 
a reasonable result with &oo = 1.27(3) [dotted curve and inset of figure 6(b)]. Given 
large uncertainties expected in the nonlinear fit of i^singha, especially for later times to, 
here we do not single out either possibility but instead provide a single final estimate 
b = 1.22(8) that covers both confidence intervals. In any case, our results indicate 
that Singha's function (9) indeed explains the global form of the temporal correlation 
function in our Eden model. This implies in particular that positive correlation would 
remain forever, i.e., lim Ct(i,to) ^ tJ > 0, which was also inferred for the circular 

i— foo 

interfaces in the liquid-crystal experiment [7]. Note that, for the flat interfaces, the 



Circular interfaces in an off-lattice Eden model 

(a) 10°^ . ^ (b) 1 



-^ 



I 



10 



10 




o ?„ = 24 



/„ = 322 



Xj = 1 -t^lt 



10 



10 




Figure 7. Short-time behavior of the temporal correlation function C't{t,to). 

(a) X2 = 1 - (^^(t, to)/(XGUE)c against xi = 1 - to/t with C[{t,tQ) = 

Ct(i, io)/(r^to*)"'^ for different reference times to- (b) Same data in the different 

2/3 

axes X2X-^ ' vs x\. The dashed lines in both panels show the right hand 

side of equation (10) with Rjl = 0.68. The dashed-dotted and dotted lines 

in the panel (b) indicate the value and the confidence interval, respectively, of 

3b2/3/25/3r(2/3) = 0.80(3). 



temporal correlation was shown to decay to zero as Ct(i, to) ~ t ^'^ for fixed to, both 
numerically [40] and experimentally [7]. Concerning the other limit t — >■ tp, equations 
(8) and (9) with Ct(io,io) = Witof ~ (rio)'/'(x'GUE>c yield 



Ct(t,io)^(rV)'/'(x')c 



1-^ 
2 



1-^ 
t 



2/3' 



(i - to < to) (10) 



with R/2 « 3fe2/3/25/3r(2/3) [39]. Our data indeed confirm this functional form 
(figure 7), except that they indicate R/2 « 0.68 rather than R/2 « 3b'^^^/2^^^T(2/S) = 
0.80(3). This implies that Singha's function (9) is not quantitatively precise in the 
short-time regime t ~ to ^ to oi the temporal correlation function. 

In relation to the liquid-crystal experiment, our data for the Eden model do 
not require the extra multiplier c(io) in equation (8), which was necessary for the 
experimental data [7]. This supports the author's speculation [7] that the need for the 
coefficient c{to) for the experiment results from microscopic dynamics of the liquid- 
crystal convection decoupled from the interface growth, which is obviously missing 
in the Eden model. With this difference in mind, the liquid-crystal experiment gave 
b w 0.8(1) and R/2 w 0.6, to be compared with b = 1.22(8) and R/2 w 0.68 in our 
Eden model. This suggests that the parameter b is not a universal quantity of the 
KPZ class. 

Next we turn our attention to another aspect of the temporal correlation, namely 
the first-passage property characterized by the temporal persistence probability 
P±{t,to). For the fluctuating interfaces, P±{t,to) is defined as the joint probability 
that the interface fiuctuation Sh{x,t) = h{x,t) — (h) at a fixed location is positive 
(negative) at time to and maintains the same sign until time t [39,40]. Here we 
distinguish the positive and negative fluctuations denoted by + and — , respectively, 
because Kallabis and Krug numerically found power-law decay P±{t,to) ^ t~^^ with 
different exponent values, 6*+ < 6*^ if A > and opposite in the other case, in the 
case of flat interfaces [40] . This asymmetry was also confirmed in the liquid-crystal 
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(a) 10' 




Figure 8. Temporal persistence probability P±(t, to) for the positive (a) and 
negative (b) fluctuations, with different to as shown in the legends. The insets 
show the running exponents — 6-|-(t, to) = d[logP-|-(t, to)]/d[log{t/to)]. The 
dashed Unes show the estimated values of the persistence exponents, 8^ = 0.81(3) 
and 9- =0.77(4). 



experiment for the flat interfaces, while no asymmetry was identified for the circular 
case [7]. Our numerically produced circular interfaces also support this claim, as 
summarized in figure 8. In the main panels, we find power-law decay 

P±(t,io)^(tAo)"'± (11) 

for sufficiently large t/tQ. This is further confirmed in the insets, where the running 
exponents ~9± {t, to) = d[log P± (t, to)]/d[log{t/to)] with different to overlap reasonably 
well when plotted against t/to and converge to constants for large abscissae. Fitting 
the data within this plateau regime (not only for the three tg in figure 8), we obtain 
9+ = 0.81(3) and 6'_ = 0.77(4). This is in good agreement with the experimentally 
found values of 6*+ = 0.81(2) and 6'_ = 0.80(2) for the circular interfaces and 
significantly different from those for the flat interfaces, 9+ = 1.35(5) and 9- = 1.85(10) 
[7]. Our result confirms in particular the absence of the asymmetry between the 
positive and negative fluctuations for the circular interfaces, i.e., 0+ = 0_ in our 
precision. The asymmetry 9^ ^ 9^ present for the flat interfaces has been attributed 
to the nonlinear term of the KPZ equation (1) [40]; understanding how this effect is 
cancelled for the circular interfaces is an interesting problem left for future studies. 

The notion of the persistence can also be considered in space [41,42]. It is then 
characterized by the spatial persistence probability Pj^ {l;t), which is the probability 
that a fluctuation remains positive or negative over length / in a spatial profile of 
the interfaces at time t. Although a few analytical and numerical studies have shown 
nontrivial character of the spatial persistence for the stationary interfaces [41,42], for 
the growing interfaces it has been studied only experimentally so far, in the liquid- 
crystal experiment [7] and in an earlier experiment on paper smoldering [43] which also 
exhibited the KPZ scaling exponents [44]. Interestingly, the two experiments showed 
different results: exponential decay for the liquid-crystal turbulence and power-law 
decay P^ {l;t) ^ Z~^/^ for the paper smoldering. The spatial persistence is therefore 
measured in our Eden model and yields the results shown in figure 9. We then find 
clear exponential decay for both positive and negative fluctuations, 

'(^'(M)~e-''i^'^, (12) 
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Figure 9. Spatial persistence probability P^ {l\t) for the positive (a) and 
negative (b) fluctuations, with different t as shown in the legends. The length 
I is shown in the rescaled unit f = {A//2)(rt)~^". The dashed lines are guides 



for the eyes showing the estimated values of k\ and k. 



Jo) 



(s) . 

with values of k±. independent of t when plotted against the dimensionless length scale 
( = (j4//2)(rt)^^/^. Measuring them at the latest two times at which the interface 
profiles were recorded, namely t = 1717 and 1999, we obtain estimates of k^ = 0.90(2) 

and K_ = 0.89(4) for our Eden model. They are in excellent agreement with the 
values numerically obtained for the temporal persistence of GUE Dyson's Brownian 
motion, k_^' = 0.90(8) and k_ = 0.90(6) [7], which is expected to be equivalent to 
the spatial profile of the curved KPZ-class interfaces. This is therefore evidence for 
universal spatial persistence in the growth regime of the KPZ class, which in fact 
concerns infinite-point correlation in the spatial profile. It is also noteworthy that the 
positive and negative fluctuations show no asymmetry in this spatial persistence either. 
Similar values were also reported in the liquid-crystal experiment, k^ = 1.07(8) and 

K_ = 0.87(6) [7], albeit with a small difference in the former which is probably due 
to finite-time effect but should be clarified by further study. In passing, for the fiat 
interfaces, the liquid-crystal experiment showed k^^' = 1.9(3) and k_' — 2.0(3) [7], 
which are significantly different from the values of the circular case. 

We finally study extreme-value statistics for the interface fiuctuations with our 
numerical data. Specifically, we focus on the maximal height -ffmax measured with 
respect to a fictitious substrate that extends from the origin [figure 10(a)], which 
is known to obey, asymptotically, the GOE Tracy- Widom distribution by studies 
of solvable models [36,45-48]§. This is the same distribution as for the one-point 
distribution of the flat interfaces and differs from the GUE Tracy- Widom distribution 
for that of the circular interfaces. For our Eden model, we measure this maximal 
height directly by ilmax = max y'^ , where i denotes the index of the active particles and 
{x',y') the coordinates defined with respect to the substrate [figure 10(a)]. Statistics 
of ffmax is improved by rotating the substrate, or the frame {x' , y'), arbitrarily around 
the origin. The obtained histograms are shown in figure 10(b) for the rescaled 
maximal height gmax = (-ffmax — Voot)/{rt)^^^, which clearly confirm the asymptotic 

§ More precisely, the conventional deflnition for the GOE Tracy- Widom random variable should be 
multiplied by 2~^" to agree with the naturally rescaled maximal height [36,45-48]. Our definition 
of the variable XGOE includes this factor 2~^'^. 
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Figure 10. Maximal height distribution, (a) Definition of the maximal height 
^max- The origin of the frame (x', y') is fixed at the position of the initial seed 
of the cluster. Shown in purple are the active particles of a cluster at t = 1999, 
which is also displayed in figure 1(b). (b) Histogram of the rescaled maximal 

height (jfnlx = (-ffmax - f oo<)/(rt)l/^ for i = 82 (blue circles) and t = 1999 
(red diamonds). The dashed and dotted lines indicate the GOE and GUE Tracy- 
Widom distributions, respectively. Note that the conventional definition of the 
GOE TW random variable is multiplied by 2~^'^, in view of the theoretical results 
for the solvable models [36,45-48]. (c) Finite-time corrections in the nth-order 
cumulants. Shown are the differences between the cumulants of the measured 
Qmax and the theoretical values for the GOE Tracy- Widom distribution, (d) 
Correction in the first-order cumulant [same data as (c)] shown in the double 
logarithmic scales. The dashed line is a guide for the eyes indicating the slope 
-2/3. 



GOE Tracy-Widom distribution for this quantity. Moreover, in figure 10(c,d), we 
measure the finite-time corrections in the nth-order cumulants ((gmax)")(. and find 
(gmax)(, — (xgoe)c ^ t^'^/'^ for the mean, in contrast with i~^/'^ found for the fiquid- 
crystal experiment [7] . Since the same set of the different exponents has been found 
for the finite-time correction in the one-point distribution, the author believes that 
the leading correction terms t^^/^ for the mean of both distributions vanish for the 
same reason in our Eden model. 
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4. Summary 

In the present paper, we have introduced an off-lattice variant of the Eden model and 
determined the statistical properties of its circular interfaces, as a numerical case study 
for the (1 + l)-diniensional KPZ universality class. Besides confirming the universal 
distribution function and the spatial correlation function in agreement with the 
analytical results on the solvable models, namely the GUE Tracy- Widom distribution 
and the Airy2 covariance, our particular focus has been put on the unsolved statistical 
properties such as the temporal correlation function and the persistence probability 
in space and time. Good agreement is then found with the recent experiment on the 
growing interfaces in the liquid-crystal turbulence, as far as characteristic quantities 
are compared such as the temporal persistence exponent 9± and the exponential 
decay rate Kj_ for the spatial persistence. In particular, unlike the flat interfaces, 
we have demonstrated that 9+ = 9^ holds for the circular case, supporting one of 
the conclusions reached in the liquid-crystal experiment. The temporal correlation 
function has also been shown to fit Singha's approximative functional form (8) 
and (9) without the small modification required for the liquid- crystal experiment. 
This indicates that the circular interfaces have long-lasting temporal correlation, i.e., 
Ct{t,to) > for large t, presumably even in the limit t — > oo unlike the flat interfaces. 
These agreements with the experimental results for the circular case and the sharp 
contrast to the flat interfaces (see table 1 of [7] for a summary) further emphasize the 
geometry-dependent universality of the (1 + l)-dimcnsional KPZ class. 

In view of these experimental and numerical results for the universality beyond 
the analytically solved quantities, it would be important to provide firmer theoretical 
grounds for such unsolved universal statistical properties, especially those related to 
the temporal correlation. This would call for several complementary approaches: 
attempts to provide analytical solutions for those quantities in solvable models, 
refinement of approximative theoretical evaluation performed by Singha [39], and 
recent developments in the application of renormalization group techniques [49- 
51] are also intriguing. Combination of such theoretical approaches as well as 
further experimental and numerical case studies would surely help understand this 
prominent out-of-equilibrium universality of the KPZ class, which governs the general 
phenomenon of the growing interfaces. 
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